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NOTATION: In the problems treated in this chapter the magnetic field B is always 
closely equal to the applied field B a , so that we write B for B a in most instances. 




Figure 1 Characteristic magnetic susceptibilities of diamagnetic and paramagnetic substances. 








CHAPTER 14 : DIAMAGNETISM AND PARAMAGNETISM 


■ 


Magnetism is inseparable from quantum mechanics, for a strictly classical 
system in thermal equilibrium can display no magnetic moment, even in a 
magnetic field. The magnetic moment of a free atom has three principal 
sources: the spin with which electrons are endowed; their orbital angular mo¬ 
mentum about the nucleus; and the change in the orbital moment induced by 
an applied magnetic field. 

The first two effects give paramagnetic contributions to the magnetization, 
and the third gives a diamagnetic contribution. In the ground Is state of the 
hydrogen atpm the orbital moment is zero, and the magnetic moment is that of 
the electron spin along with a small induced diamagnetic moment. In the Is 2 
state of helium the spin and orbital moments are both zero, and there is only an 
induced moment. Atoms with filled electron shells have zero spin and zero 
orbital moment: these moments are associated with unfilled shells. 

The magnetization M is defined as the magnetic moment per unit volume. 
The magnetic susceptibility per unit volume is defined as 


(cgs) x = Y ; 


(si) 


X ~ 


J*qM 

B 


( 1 ) 


where B is the macroscopic magnetic field intensity. In both systems of units x 
is dimensionless. We shall sometimes for convenience refer to M/B as the sus¬ 
ceptibility without specifying the system of units. 

Quite frequently a susceptibility is defined referred to unit mass or to a 
mole of the substance. The molar susceptibility is written as the magnetic 
moment per gram is sometimes written as cr. Substances with a negative mag¬ 
netic susceptibility are called diamagnetic. Substances with a positive suscepti¬ 
bility are called paramagnetic, as in Fig. 1. 

Ordered arrays of magnetic moments are discussed in Chapter 15; the 
arrays may be ferromagnetic, ferrimagnetic, antiferromagnetic, helical, or 
more complex in form. Nuclear magnetic moments give rise to nuclear 
paramagnetism. Magnetic moments of nuclei are of the order of 10 -3 times 
smaller than the magnetic moment of the electron. 


LANGEVIN DIAMAGNETISM EQUATION 

Diamagnetism is associated with the tendency of electrical charges par¬ 
tially to shield the interior of a body from an applied magnetic field. In electro¬ 
magnetism we are familiar with Lenz’s law: when the flux through an electrical 
circuit is changed, an induced current is set up in such a direction as to oppose 
the flux change. 
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In a superconductor or in an electron orbit within an atom, the induced 
current persists as long as the field is present. The magnetic field of the induced 
current is opposite to the applied field, and the magnetic moment associated 
with the current is a diamagnetic moment. Even in a normal metal there is a 
diamagnetic contribution from the conduction electrons, and this diamag¬ 
netism is not destroyed by collisions of the electrons. 

The usual treatment of the diamagnetism of atoms and ions employs the 
Larmor theorem: in a magnetic field the motion of the electrons around a 
central nucleus is, to the first order in B, the same as a possible motion in the 
absence of B except for the superposition of a precession of the electrons with 
angular frequency 


(CGS) a) = eB/2mc ; 


(SI) a) = eB!2m . 


( 2 ) 


If the field is applied slowly, the motion in the rotating reference system will be 
the same as the original motion in the rest system before the application of the 
field. 

If the average electron current around the nucleus is zero initially, the 
application of the magnetic field will cause a finite current around the nu¬ 
cleus. The current is equivalent to a magnetic moment opposite to the applied 
field. It is assumed that the Larmor frequency (2) is much lower than the fre¬ 
quency of the original motion in the central field. This condition is not satisfied 
in free carrier cyclotron resonance, and the cyclotron frequency is twice the 
frequency (2). 

The Larmor precession of Z electrons is equivalent to an electric current 


(SI) 


= (charge)(revolutions per unit time) = (— Ze)(— 

'27 



eB \ 

Utt 

2m / 


( 3 ) 


The magnetic moment /x of a current loop is given by the product 
(current) X (area of the loop). The area of the loop of radius p is 7rp 2 . We have 


(SI) 




Z e 2 B 
4 m 


<P 2 > ; 


(CGS) n = — 


Z£B 

4mc 2 


<p 2 > • (4) 


Here (p 2 ) = (x 2 ) + (y 2 ) is the mean square of the perpendicular distance of the 
electron from the field axis through the nucleus. The mean square distance of 
the electrons from the nucleus is (r 2 ) = (x 2 ) + (y 2 ) + (z 2 ). For a spherically 
symmetrical distribution of charge we have (x 2 ) = (y 2 ) = (z 2 ), so that (r 2 ) = 

Kp 2 >. 

From (4) the diamagnetic susceptibility per unit volume is, if N is the 
number of atoms per unit volume, 


X = 


Nfi 


NZe 2 
6 me 2 


<r 2 > ; 


(CGS) 


B 


(5) 
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(SI) 

B 6m 

This is the classical Langevin result. 

The problem of calculating the diamagnetic susceptibility of an isolated 
atom is reduced to the calculation of (r 2 ) for the electron distribution within the 
atom. The distribution can be calculated by quantum mechanics. 

Experimental values for neutral atoms are most easily obtained for the 
inert gases. Typical experimental values of the molar susceptibilities are the 
following: 


He Ne Ar Kr Xe 

Xm in CGS in 10 -6 cm 3 /mole: —1.9 -7.2 —19.4 -28.0 -43.0 


■M 




iHHHH 


In dielectric solids the diamagnetic contribution of the ion cores is de¬ 
scribed roughly by the Langevin result. The contribution of conduction elec¬ 
trons is more complicated, as is evident from the de Haas-van Alphen effect 
discussed in Chapter 9. 


QUANTUM THEORY OF DIAMAGNETISM OF MONONUCLEAR SYSTEMS 


From (G. 18) the effect of a magnetic field is to add to the hamiltonian the 
terms 


ip.fi P 2 

X =-(V- A + A-V) + — T 

2mc 2 me 2 


( 6 ) 


for an atomic electron these terms may usually be treated as a small perturba¬ 
tion. If the magnetic field is uniform and in the z direction, we may write 


A x = -it/B , A y = ixB , A z = 0 , 
and (6) becomes 


iehB f d d \ 
2 me \ dy ^ dx / 


+ 


e 2 B 2 
8 me 2 


(x 2 + y 2 ) 


(7) 

( 8 ) 


The first term on the right is proportional to the orbital angular momen¬ 
tum component L z if r is measured from the nucleus. In mononuclear systems 
this term gives rise only to paramagnetism. The second term gives for a spheri¬ 
cally symmetric system a contribution 














4m 


by first-order perturbation theory. The associated magnetic moment is diamag¬ 
netic: 


ar = e 2 (r 2 ) 

BB 6 me 2 


( 10 ) 


in agreement with the classical result. 


PARAMAGNETISM 

Electronic paramagnetism (positive contribution to %) * s found in: 

L Atoms, molecules, and lattice defects possessing an odd number of 
electrons, as here the total spin of the system cannot be zero. Examples: free 
sodium atoms; gaseous nitric oxide (NO); organic free radicals such as 
triphenylmethyl, C(C 6 H 5 ) 3 ; F centers in alkali halides, 

2. Free atoms and ions with a partly filled inner shell; transition elements; 
ions isoelectronic with transition elements; rare earth and actinide elements. 
Examples: Mn 2+ , Gd 3+ , U 4+ . Paramagnetism is exhibited by many of these 
ions even when incorporated into solids, but not invariably. 

■3. A few compounds with an even number of electrons, including molecu¬ 
lar oxygen and organic biradicals. 

4. Metals. 


QUANTUM THEORY OF PARAMAGNETISM 


The magnetic moment of an atom or ion in free space is given by 




(ii) 


where the total angular momentum &J is the sum: of the: orbital hh and spin hS 
angular momenta. 

The constant y is the ratio of the magnetic moment to the angular momen¬ 
tum; y is called the gyromagnetie ratio or magnetogyric ratio. For electronic 
systems a quantity g called the g factor or the spectroscopic, splitting factor is 
defined by 

g^s = - yk ■ (12) 

For an electron spin g = 2.0023, usually taken as 2.00, For a free atom the g 
factor is given by the Eande equation 

/(/ + 1) + S(S + 1) - L(L + 1) 

8 = 1 + 


mi + 1 ) -- 


(13) 




14 Diamagnetism and Paramagnetism 


421 



2fiB 


m s 

2 


Hz 

~H 


2 


H 


Figure 2 Energy level splitting for one electron 
in a magnetic field B directed along the positive z 
axis. For an electron the magnetic moment fx is 
opposite in sign to the spin S, so that /x = 
-gfx B S. In the low energy state the magnetic 
moment is parallel to the magnetic field. 



Figure 3 Fractional populations of a two-level 
system in thermal equilibrium at temperature T 
in a magnetic field B. The magnetic moment is 
proportional to the difference between the two 
curves. 


The Bohr magneton fig is defined as eh/2mc in CGS and ehl2m in SI. It is 
closely equal to the spin magnetic moment of a free electron. 

The energy levels of the system in a magnetic field are 

U = -M ’ B = , (14) 

where nij is the azimuthal quantum number and has the values /, J — 1, . . . , 
—J. For a single spin with no orbital moment we have mj = ±i and g = 2, 
whence U = ±/jl b B. This splitting is shown in Fig. 2. 

If a system has only two levels the equilibrium populations are, with 

r = k B T 9 

exp (fxB/r) 


Ah 

N 

^ 2 _ 

N 


exp(/jiB/T) + exp (—fiB/r) 
exp(— (iB/r) 


(15) 


(16) 


exp (/jlB/t) + exp (—/xB/r) 

here N 2 are the populations of the lower and upper levels, and 
N = N\ + N 2 is the total number of atoms. The fractional populations are plot¬ 
ted in Fig. 3. 

The projection of the magnetic moment of the upper state along the field 
direction is —fi and of the lower state is fi. The resultant magnetization for N 
atoms per unit volume is, with x = ixB/k B T , 


M = (Ni ~ N 2 )jjl = Nfx 


e v — e 
e x + e - 


= N/jl tanh x . 


(17) 


For x < 1, tanh x = x, and we have 

M = Nfx{jxBlk B T) . (18) 

In a magnetic field an atom with angular momentum quantum number } 
has 2/ + 1 equally spaced energy levels. The magnetization (Fig. 4) is given by 

M = NgJfi B Bj(x) , (x - gJn B B/k B T) , (19) 
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Figure 4 Plot of magnetic moment versus BIT for spherical samples of (I) potassium chromium 
alum, (II) ferric ammonium alum, and (III) gadolinium sulfate octahydrate. Over 99.5% magnetic 
saturation is achieved at 1.3 K and about 50,000 gauss. (5 T). After W. E. Henry. 


where the Brillouin function Bj is defined by 


2 / + 1 

Bj(x) = ——— ctnh 


(2 1 ± D* 

2 J 


-ctnh — 


2 / 


2 / 


Equation (17) is a special case of (20) for J = h 
For x 1, we have 


ctnh x = — 
x 


45 


+ 


and the susceptibility is 

M _ NJ(J + l)gVI AfoVI C 
B 3 k B T 3k B T T 

Here p is the effective number of Bohr magnetons, defined as 


( 20 ) 


( 21 ) 


( 22 ) 


v * gUU + D ] 1/2 


( 23 ) 
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Figure 5 Plot of 1/* os T for a gadolinium salt, Gd(C 2 H 5 S0 4 )3 • 9H 2 0, The straight line is the 
Curie law. (After L, C, Jackson and B. Kamerlingh Onnes.) 


The constant C is known as the Curie constant. The form (19) is known as the 
Curie-Brillouin law, and (22) is known as the Curie law, Results for the para¬ 
magnetic ions in a gadolinium salt are shown in'.Fig, 5, 

Mare Earth Ions 

The ions of the rare-earth elements have closely similar chemical proper¬ 
ties, and their chemical separation in tolerably pure form was accomplished 
only long after their discovery. Their magnetic properties are fascinating: the 
ions exhibit a systematic variety and intelligible complexity. The chemical 
properties of the trivalent ions are similar because the outermost electron shells 
are identically in the Ss^Bp 6 configuration, like neutral xenon. In lanthanum, 
just before the rare earth group begins, the 4/ shell is empty; at cerium there is 
one 4 f electron, and the number of 4/electrons increases steadily through the 
group until we have 4/ 13 at ytterbium and the filled shell if 14 at lutecium. The 
radii of the trivalent ions contract fairly smoothly as we go through the group 
from 1,11 A at cerium to 0,94 A at ytterbium. This is the famous ‘lanthanide 
contraction.” What distinguishes the magnetic behavior of one ion species from 
another is the number of 4/electrons compacted in the inner shell with a radius 
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of perhaps 0.3 A. Even in the metals the 4/ core retains its integrity and 
its atomic properties: no other group of elements in the periodic table is as 
interesting. 

The preceding discussion of paramagnetism applies to atoms that have a 
(2/ + I)~fold degenerate ground state, the degeneracy being lifted by a mag¬ 
netic field. The influence of all higher energy states of the system is neglected. 
These' assumptions -appear to be satisfied by a number of rare-earth ions, 
Table 1. The calculated magneton numbers are obtained with g values from the 
Lande result (13) and the ground-state level assignment predicted by the Hund 
theory of spectral terms. The discrepancy between the experimental magneton 
numbers and those calculated on these assumptions is quite marked for Eu 3+ 
and Sm 3+ ions. For these ions it is necessary to consider the influence of the 
high states of the L-S multiple^ as the intervals between successive states of 
the multiple! are not large compared to hgT at room temperature. A multiplet 
is the set of levels of different J values arising out of a given L and S, The levels 
of a multiple! are split by the spin-orbit interaction. 

Hund Rules 

The Hund rules as applied to electrons in a given shell of hn atom affirm 
that electrons will occupy orbitals in such a way that the ground state is charac¬ 
terized by the following: 

;L The maximum value of the total spin S allowed by the exclusion 
principle; 

2; The maximum value of the orbital angular momentum L consistent with 
this value: of S; 

3* The value of the total angular momentum / is equal to |L - S j when the 
shell is less than half full and to L + 5 : when the: shell is more than half full. 
When the shell is just half full, the application of the first rule gives L — 0, so 
that J = S. 

The first Hund rule has its origin In the exclusion principle and the cou¬ 
lomb repulsion between electrons. The exclusion principle prevents two elec¬ 
trons of the same spin from being at the same place at the same time. Thus 
electrons of the same spin are kept apart, further apart than electrons of oppo¬ 
site spin. Because of the coulomb interaction the energy of electrons of the 
same spin is lower—the average potential energy is less positive for parallel 
spin than for antiparallel spin. A good example is the ion Mn 2+ . This ion has 
five electrons in the 3 d shell, which is therefore half-filled. The spins can all be 
parallel if each electron enters a different orbital, and there are exactly five 
different orbitals available, characterized by the orbital quantum numbers 
m L = 2, 1, (X - I,. -2, These will be occupied by one electron each. We expect 
S = I, and because = 0 the only possible value of L is 0, as observed. 
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Table 1 Effective magneton numbers p for trivalent lanthanide group ions 
(Near room temperature) 


Ion 

Configuration 

Basic level 

p(calc) = 

g UU + i)] 172 

p(exp), 

approximate 

Ce 3+ 

4f5s 2 p 6 

2 F5/2 

2.54 

2.4 

p r 3+ 

4/ 2 5.s 2 p 6 

3 h 4 

3.58 

3.5 

Nd 3+ 

4f 3 5s 2 p 6 

4 h/2 

3.62 

3.5 

Pm 3+ 

4/ 4 5s 2 p 6 

% 

2.68 

— 

Sm 3+ 

4/ 5 5s 2 p 6 

6 Hs/2 

0.84 

1.5 

Eu 3+ 

4/ 6 5s 2 p 6 

7 Fo 

0 

3.4 

Gd 3+ 

4f5s 2 p 6 

8c 

07/2 

7.94 

8.0 

Tb 3+ 

4f5 s 2 p 6 

7 F 6 

9.72 

9.5 

Dy 3+ 

4f 9 5sY 

6 Hl5/2 

10.63 

10.6 

Ho 3+ 

4f°5sY 

5t 

^8 

10.60 

10.4 

Er 3+ 

4f n 5sY 

4t 

115/2 

9.59 

9.5 

Tm 3+ 

4f 12 5sY 

3 h 6 

7.57 

7.3 

Yb 3+ 

4f 13 5sY 

2 F 7 /2 

4.54 

4.5 


The second Hund rule is best approached by model calculations. Pauling 
and Wilson, 1 for example, give a calculation of the spectral terms that arise from 
the configuration p 2 . The third Hund rule is a consequence of the sign of the 
spin-orbit interaction; For a single electron the energy is lowest when the spin 
is antiparallel to the orbital angular momentum. But the low energy pairs m Ly 
m$ are progressively used up as we add electrons to the shell; by the exclusion 
principle when the shell is more than half full the state of lowest energy neces¬ 
sarily has the spin parallel to the orbit. 

Consider two examples of the Hund rules; The ion Ce 3+ has a single f 
electron; an f electron has l = 3 and s = f. Because the/shell is less than half 
full, the J value by the preceding rule is \L — S| = L — % = f. The ion Pr 3+ has 
two /electrons: one of the rules tells us that the spins add to give S = 1. Both/ 
electrons cannot have mi = 3 without violating the Pauli exclusion principle, so 
that the maximum L consistent with the Pauli principle is not 6, but 5. The J 
value is |L — S| = 5 — 1 = 4. 

Iron Group Ions 

Table 2 shows that the experimental magneton numbers for salts of the iron 
transition group of the periodic table are in poor agreement with (18). The 
values often agree quite well with magneton numbers p = 2[S(S + 1)] 1/2 calcu- 

J L. Pauling and E. B. Wilson, Introduction to quantum mechanics , McGraw-Hill, 1935, 
pp. 239-246. 
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Table 2 Effective magneton numbers for iron group ions 


Ion 

Config¬ 

uration 

Basic 

level 

p(calc) = 

glJU + l)] 1 ' 2 

p(calc) = 

2[S(S + l)] 1 ' 2 

p(exp) s 


Ti 3+ , V 4+ 

3 d l 

2 C>3/2 

1.55 

1.73 

1.8 

y3+ 

3 d 2 

3 f 2 

1.63 

2.83 

2.8 

Cr 3+ , V 2+ 

3 d 3 

4 Fm 

0.77 

3.87 

3.8 

Mn 3+ , Cr 2+ 

3d 4 

5 D 0 

0 

4.90 

4.9 

Fe 3+ , Mn 2+ 

3d 5 

6 Ss/ 2 

5.92 

5.92 

5.9 

Fe 2+ 

3d 6 

5 d 4 

6.70 

4.90 

5.4 

Co 2+ 

3d 7 

4 f 9/2 

6.63 

3.87 

4.8 

Ni 2+ 

3d 8 

3 f 4 

5.59 

2.83 

3.2 

Cu 2+ 

3d 9 

2 F>5/2 

3.55 

1.73 

1.9 


a Representative values. 


lated as if the orbital moment were not there at all. We say that the orbital 
moments are quenched. 

Crystal Field Splitting 

The difference in behavior of the rare earth and the iron group salts is that 
the 4 f shell responsible for paramagnetism in the rare earth ions lies deep 
inside the ions, within the 5 s and 5 p shells, whereas in the iron group ions the 
3 d shell responsible for paramagnetism is the outermost shell. The 3 d shell 
experiences the intense inhomogeneous electric field produced by neighboring 
ions. This inhomogeneous electric field is called the crystal field. The interac¬ 
tion of the paramagnetic ions with the crystal field has two major effects: the 
coupling of L and S vectors is largely broken up, so that the states are no longer 
specified by their J values; further, the 2L + 1 sublevels belonging to a given L 
which are degenerate in the free ion may now be split by the crystal field, as in 
Fig. 6. This splitting diminishes the contribution of the orbital motion to the 
magnetic moment. 

Quenching of the Orbital Angular Momentum 

In an electric field directed toward a fixed nucleus, the plane of a classical 
orbit is fixed in space, so that all the orbital angular momentum components L r , 
L y7 L z are constant. In quantum theory one angular momentum component, 
usually taken as L z , and the square of the total orbital angular momentum L 2 are 
constant in a central field. In a noncentral field the plane of the orbit will move 
about; the angular momentum components are no longer constant and may 
average to zero. In a crystal L z will no longer be a constant of the motion, 
although to a good approximation L 2 may continue to be constant. When L z 
averages to zero, the orbital angular momentum is said to be quenched. The 
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Figure 6 Consider an atom with orbital angular momentum L = 1 placed in the uniaxial crystalline 
electric field of the two positive ions along the z axis. In the free atom the states m h = ±1, 0 have 
identical energies—they are degenerate. In the crystal the atom has a lower energy when the 
electron cloud is close to positive ions as in (a) than when it is oriented midway between them, as 
in (b) and (c). The wavefunctions that give rise to these charge densities are of the form zf(r), xf(r) 
and yf(r) and are called the p z> p z , p y orbitals, respectively. In an axially symmetric field, as shown, 
the p x and p y orbitals are degenerate. The energy levels referred to the free atom (dotted line) are 
shown in (d). If the electric field does not have axial symmetry, all three states will have different 
energies. 


magnetic moment of a state is given by the average value of the magnetic 
moment operator ^t B (L + 2S). In a magnetic field along the z direction the 
orbital contribution to the magnetic moment is proportional to the quantum 
expectation value of L z ; the orbital magnetic moment is quenched if the me¬ 
chanical moment L z is quenched. 

When the spin-orbit interaction energy is introduced, the spin may drag 
some orbital moment along with it. If the sign of the interaction favors parallel 
orientation of the spin and orbital magnetic moments, the total magnetic mo¬ 
ment will be larger than for the spin alone, and the g value will be larger than 2. 
The experimental results are in agreement with the known variation of sign of 
the spin-orbit interaction: g > 2 when the 3 d shell is more than half full, g = 2 
when the shell is half full, and g < 2 when the shell is less than half full. 

We consider a single electron with orbital quantum number L = 1 moving 
about a nucleus, the whole being placed in an inhomogeneous crystalline elec¬ 
tric field. We omit electron spin. 

In a crystal of orthorhombic symmetry the charges on neighboring ions 
will produce an electrostatic potential (p about the nucleus of the form 

e<p = Ax 2 + By 2 - (A + B)z 2 , (24) 

where A and B are constants. This expression is the lowest degree polynomial 
in x , y , z which is a solution of the Laplace equation V 2 (p = 0 and compatible 
with the symmetry of the crystal. 
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In free space the ground state is three-fold degenerate, with magnetic 
quantum numbers m L ~ 1, 0, — L In a magnetic field these levels are split by 
energies proportional to the field B, and it is this field-proportional splitting 
which is responsible for the normal paramagnetic susceptibility of the ion. In 
the crystal the picture may be different, We take as the three wavefunotions 
associated with the unperturbed ground state of the ion 

U x = xf(r) ; U y - yf(r) ; U z - zf(r) . (25) 

These wavefunetions are orthogonal, and we assume that they are normalized. 
Each of the Its mn be shown to have the property 

^Ui - UL + Wi = 2 Vi , (26) 

where i£ 2 is the operator for the square of the orbital angular momentum, in 
units of h. The result (26) confirms that the selected wavefunetions are in fact p 
functions, having L — 1. 

We observe now that the U*$ are diagonal with respect to the perturbation, 
as by symmetry the nondiagonal elements vanish: 

<14M14> = (U x \e<p\U z ) = (V y \e<p\U z ) - 0 . (27) 

Consider for example 

(U x \ef\Uy) — / xy\f(r)\ 2 {Ax 2 T By 2 ~~ (A- + dx dy dz ; (28) 

the integrand is an odd function of x (and also of y) and therefore the integral 
must be zero. The energy levels are then given by the diagonal matrix 
elements: 


(U x \e<p\U x ) - f [/(r)| 2 {A# + Bifx 2 - (A + B)z 2 xF} dx dy dz 

= Mh - h) , 


(29) 


where 

h = / |/(r)| 2 x 4 dx dy dz ; 1 2 = / \f(r)\ 2 xFy 2 dx dy dz . 

In addition, 

(U y \e<p\U y ) - B (l| - h) ; (U z \e<p\U z ) - —(A + B)(h - h) . 

Hie three eigenstates in the crystal field are p functions with their angular lobes 
directed along each of the x, f, z axes,, respectively. 

The orbital moment of each of tie level's is zero, because 


(VMM - WiL z \U V ) = (U Z \L Z \U Z ) = 0 . 

This effect is known as quenching. The level still has a definite total angular 
momentum, since it 2 is diagonal and gives L — : 1, but the spatial components 
of the angular momentum are not constants of the motion and their time aver¬ 
age is zero in the first approximation. Therefore the components of the orbital 
magnetic moment also vanish in the same approximation. The role of the crystal 
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field in the quenching process is to split the originally degenerate levels into 
nonmagnetic levels separated by energies §> pH, so that the magnetic field is a 
small perturbation in comparison with the crystal field. 

At a lattice site of cubic symmetry there is no term in the potential of the 
form (24), that is, quadratic in the electron coordinates. Now the ground state 
of an ion with one p electron (or with one hole in a p shell) will be triply 
degenerate. However, the energy of the ion will be lowered if the ion displaces 
itself with respect to the surroundings, thereby creating a noncubic potential 
such as (24). Such a spontaneous displacement is known as a Jahn-Teller effect 
and is often large and important, particularly with the Mn 3+ and Cu 2+ ions 2 
and with holes in alkali and silver halides. 

Spectroscopic Splitting Factor 

We suppose for convenience that the crystal field constants, A, B are such 
that U x ” xf(r) is the orbital wave function of the ground state of the atom in 
the crystal. For a spin S = i there are two possible spin states S % = ±i repre¬ 
sented by the spin functions or, /3, which in the absence of a magnetic field are 
degenerate in the zeroth approximation. The problem is to take into account 
the spin-orbit interaction energy AL* S, 

If the ground state function is = U-^a = xf(r)a in the zeroth approxima¬ 
tion, then in the first approximation, considering the AL * S interaction by 
standard perturbation theory, we have 

if, = [14 - i(k/2&i)U y ]a - HA/m^US , (30) 

where Ai is the energy difference between the U x and U y states, and A 2 is the 
difference between the U x and U z states. The term in 3 actually has only a 
second order effect on the result and may be discarded. The expectation value 
of the orbital angular momentum to the first order is given directly by 

(*W«0 = -A/A! , 

and the magnetic moment of the state as measured in the % direction is 
/%($£* + 2$j0) = |.-(A/A|) + 13$#. , 

As the separation between the levels S* = ±| in a field W is 
AE « gp s li « 2[1 - , 

the g value or spectroscopic splitting factor (12) in the z direction is 

g = 2[1 - (A/A,)] . (31) 


2 See L. Orgel, Introduction to transition metal dhemistmf:, 2nd : ecL :> Wiley, 1966; extensive 
references are given by M. D. Sturge, Phys. Rev. 149, A880 (1965), 
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Van Vleck Temperature-Independent Paramagnetism 

We consider an atomic or molecular system which has no magnetic mo¬ 
ment in the ground state, by which we mean that the diagonal matrix element 
of the magnetic moment operator p, z is zero. 

Suppose that there is a nondiagonal matrix element (s|//JO) of the magnetic 
moment operator, connecting the ground state 0 with the excited state s of 
energy A = E s - E 0 above the ground state. Then by standard perturbation 
theory the wavefunction of the ground state in a weak field (/ x z B < A) becomes 

</'& = </'o + (B/&)(s\p, z \0)ip s , (32) 

and the wavefunction of the excited state becomes 

# = “ (B/A)<0|/tJ*>*> . (33) 

The perturbed ground state now has a moment 

<0'|m*|0'> = 2B|(s|/x 3 |0)| 2 /A , (34) 

and the upper state has a moment 

<*Vsk>= -2B|( 5 |^|0)| 2 /A . (35) 

There are two interesting cases to consider: 

Case (a). A k B T. The surplus population in the ground state over the 
excited state is approximately equal to NA/2/c fi T, so that the resultant magneti¬ 
zation is 


M -, ,36) 

which gives for the susceptibility 

X = NM/Lt 2 |0>| 2 /fc e r . (37) 

Here N is the number of molecules per unit volume. This contribution is of the 
usual Curie form, although the mechanism of magnetization here is by polariza¬ 
tion of the states of the system, whereas with free spins the mechanism of 
magnetization is the redistribution of ions among the spin states. We note that 
the splitting A does not enter in (37). 

Case ( b ). A > k B T. Here the population is nearly all in the ground state, so 

that 


M = 


2NBl<s| / Lt z |0)| 2 

A 


(38) 


The susceptibility is 


2A/|(s|Mz|0)| 2 


A 


(39) 
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independent of temperature. This type of contribution it known as Van Vleck 
paramagnetism. 


COOLING BY ISENTMOPIC DEMAGNETIZATION 

The first method for attaining temperatures much below 1 K was that of 
isentropie, or adiabatic, demagnetization, 3 By its use temperatures of 10" 3 K 
and lower have been reached. The method rests on the fact that at a fixed 
temperature the entropy of a system of magnetic moments is lowered by the 
application of a magnetic field. 

The entropy is a measure of the disorder of a system: the greater the 
disorder, the higher is the entropy. In the magnetic field the moments will be 
partly lined up (partly ordered), so that the entropy is lowered by the field. The 
entropy is also lowered if the temperature is lowered, as more of the moments 
line up. 

If the magnetic:'.field can then be removed without changing the entropy of 
the spin system, the order of the spin system will look like a lower temperature 
than the same degree of order in the presence of the field. When the specimen 
is demagnetized at constant entropy, entropy can flow into the spin system only 
from the system of lattice vibrations, as in Fig, 7, At the temperatures of inter¬ 
est the entropy of the lattice vibrations is usually negligible; thus the entropy of 
the spin system will be essentially constant during adiabatic demagnetization of 
the specimen. Magnetic cooling is a one-shot operation, not cyclic. 

We first find an expression for the spin entropy of a system of N ions, each 
of spin S, at a temperature sufficiently high that the spin system is entirely 
disordered. That is, T is supposed to be much higher than some temperature A 
which characterizes the energy of the interactions (E int = k B A) tending to ori¬ 
ent the spins preferentially. Some of these interactions are discussed in Chap¬ 
ter 16. The definition of the entropy cr of a system of G accessible states is 
<t = k B In G. At a temperature so high that all of the 2S + 1 states of each ion 
are nearly equally populated, G is the number of ways of arranging N spins in 
2S + 1 states. Thus G = (2S + if, whence the spin entropy cr s is: 

crs - k B In (IS + if - Nkg In (2S T 1) . (40) 

This spin entropy is reduced by a magnetic field if the lower levels gain in 
population '.when the field; separates the $$' + I states m energy. 


3 The method was suggested by P. Debye, Ann. Physik 81, 1154 (1926); and W. F, Giauque, 
]. Am. Chem, Soc. 49, 1864 (1927). For many purposes the method has been supplanted by the 
He 3 -He 4 dilution refrigerator which operates in a continuous cycle. The He 3 atoms in solution in 
liquid He 4 play the role of atoms in a gas, and cooling is effected by “vaporization” of He 3 ;: see TP, 
Chapter 12. 
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Figure 7 During adiabatic demagnetization the total entropy of the specimen is constant. For 
effective cooling the initial entropy of the lattice should be small in comparison with the entropy of 
the spin system. 


The steps carried out in the cooling process are shown in Fig. 8. The field 
is applied at temperature T\ with the specimen in good thermal contact with 
the surroundings, giving the isothermal path ab. The specimen is then insu- 
lated (Act = 0) and the field removed; the specimen follows the constant en¬ 
tropy path be , ending up at temperature T 2 . The thermal contact at T\ is pro¬ 
vided by helium gas, and the thermal contact is broken by removing the gas 
with a pump. 


Nuclear Demagnetization 


The population of a magnetic sublevel is a function only of (iB/k B T, hence 
of BIT. The spin-system entropy is a function only of the population distribu¬ 
tion; hence the spin entropy is a function only of BIT. If is the effective field 
that corresponds to the local interactions, the final temperature T 2 reached in 
an adiabatic demagnetization experiment is 


T 2 = T\{BJB) , 


(41) 


where B is the initial field and T± the initial temperature. 

Because nuclear magnetic moments are weak, nuclear magnetic interac¬ 
tions are much weaker than similar electronic interactions. We expect to reach 
a temperature 100 times lower with a nuclear paramagnet than with an electron 
paramagnet. The initial temperature 7\ of the nuclear stage in a nuclear spin¬ 
cooling experiment must be lower than in an electron spin-cooling experiment. 
If we start at B = 50 kG and T x = 0.01 K, then i±Blk B Ti ~ 0.5, and the en- 
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7, mK 

' igure 8 Entropy for a spin £ system as a function of temperature, assuming an internal random 
magnetic field B A of 100 gauss. The specimen is magnetized isothermally along ab , and is then 
insulated thermally. The external magnetic field is turned off along be. In order to keep the figure 
on a reasonable scale the initial temperature Ti is lower than would be used in practice, and so is 
the external magnetic field. 

tropy decrease on magnetization is over 10 percent of the maximum spin en¬ 
tropy. This is sufficient to overwhelm the lattice and from (41) we estimate a 
final temperature T 2 ~ 10~ 7 K. The first 4 nuclear cooling experiment was car¬ 
ried out on Cu nuclei in the metal, starting from a first stage at about 0.02 K 
as attained by electronic cooling. The lowest temperature reached was 
1.2 X 10" 6 K. 

The results in Fig. 9 fit a line of the form of (41): T 2 = 7^(3.1/B) with B in 
gauss, so that B A = 3.1 gauss. This is the effective interaction field of the mag¬ 
netic moments of the Cu nuclei. The motivation for using nuclei in a metal is 
that conduction electrons help ensure rapid thermal contact of lattice and nu¬ 
clei at tine temperature of the first stage. The present record 5 for a spin temper¬ 
ature is 280 pK, in rhodium. 


PARAMAGNETIC SUSCEPTIBILITY OF CONDUCTION ELECTRONS 

We are going to try to show how on the basis of these statistics the fact that many 
metals are diamagnetic, or only weakly paramagnetic, can be brought into agree¬ 
ment with the existence of a magnetic moment of the electrons. 

W. Pauli, 1927 

Classical free electron theory gives an unsatisfactory account of the para¬ 
magnetic susceptibility of the conduction electrons. An electron has associated 
with it a magnetic moment of one Bohr magneton, fx B . One might expect that 

4 N. Kurti, F. N. H. Robinson, F. E. Simon, and D. A. Spohr, Nature 178, 450 (1956); for 
reviews see N. Kurti, Cryogenics 1, 2 (1960); Adv. in Cryogenic Engineering 8, 1 (1963). 

5 P. J. Hakonen et al., Phys. Rev . Lett . 70, 2818 (1993). 
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Initial magnetic field in kG 



Figure 9 Nuclear demagnetizations of copper nuclei in the metal, starting from 0.012 K and 
various fields. (After M. V. Hobden and N. Kurti.) 


the conduction electrons would make a Curie-type paramagnetic contribution 
(22) to the magnetization of the metal: M = N/JL B B/k B T. Instead it is observed 
that the magnetization of most normal nonferromagnetic metals is independent 
of temperature. 

Pauli showed that the application of the Fermi-Dirac distribution (Chap¬ 
ter 6) would correct the theory as required. We first give a qualitative explana¬ 
tion of the situation. The result (18) tells us that the probability an atom will be 
lined up parallel to the field B exceeds the probability of the antiparallel orien¬ 
tation by roughly fxBlk B T. For N atoms per unit volume, this gives a net mag¬ 
netization ^Nfj^B/ksT, the standard result. 

Most conduction electrons in a metal, however, have no possibility of 
turning over when a field is applied, because most orbitals in the Fermi sea 
with parallel spin are already occupied. Only the electrons within a range k B T 
of the top of the Fermi distribution have a chance to turn over in the field; thus 
only the fraction TIT F of the total number of electrons contribute to the suscep¬ 
tibility. Hence 


N/jl 2 B T _ Ny 
k B T T f k B Tp 


which is independent of temperature and of the observed order of magnitude. 

We now calculate the expression for the paramagnetic susceptibility of a 
free electron gas at T < Tp. We follow the method of calculation suggested by 
Fig. 10. An alternate derivation is the subject of Problem 5. 
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Total energy, kinetic + 
magnetic, of electrons 



Density of 
orbitals 


Figure 10 Pauli paramagnetism at absolute zero; the orbitals in the shaded regions in (a) are 
occupied. The numbers of electrons in the “up” and “down” band will adjust to make the energies 
equal at the Fermi level. The chemical potential (Fermi level) of the moment up electrons is equal 
to that of the moment down electrons. In (b) we show the excess of moment up electrons in the 
magnetic field. 


The concentration of electrons with magnetic moments parallel to the 
magnetic field is 


N+ = d- [ de D(e + i*B) = -1- [ de D(e) + L /xB D(e F ) , 

* J ~nB * J Q ^ 

written for absolute zero. Here \D(e + /jlB) is the density of orbitals of one 

spin orientation, with allowance for the downward shift of energy by — fiB. 
The approximation is written for k B T < e F . 

The concentration of electrons with magnetic moments antiparallel to the 
magnetic field is 


N- = — de D(e — (jlB) = — 


HB 


de D(e) —— ^xB D(e F ) . 


The magnetization is given by M = fi(N+ — N_), so that 

o 3 Nul 2 

M = M 2 D(e F )B = —!B , (42) 

ZKbI F 

with D(e F ) ~ 3N/2e F = 3N/2k B T F from Chapter 6. The result (42) gives the 
Pauli spin magnetization of the conduction electrons, for k B T < e F . 

In deriving the paramagnetic susceptibility, we have supposed that the 
spatial motion of the electrons is not affected by the magnetic field. But the 
wavefunctions are modified by the magnetic field; Landau has shown that for 
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free electrons this causes a diamagnetic moment equal to — 3 of the paramag¬ 
netic moment. Thus the total magnetization of a free electron gas is 


M 


k B T F 


B . 


(43) 


Before comparing (43) with the experiment we must take account of the 
diamagnetism of the ionic cores, of band effects, and of electron-electron inter¬ 
actions. In sodium the interaction effects increase the spin susceptibility by 
perhaps 75- -percent. 

The magnetic susceptibility is considerably higher for most transition met¬ 
als (with unfilled inner electron shells) than for the alkali metals (Fig. 11). The 
high values suggest that the density of orbitals is unusually high for transition 
metals, in agreement with measurements of the electronic heat capacity. We 
saw fii Chapter 9 how this arises from band theory. 


SUMMARY 
(In CGS Units) 

* The diamagnetic susceptibility of N atoms of atomic number Z is x = 
— Ze 2 N(r 2 )/6mc 2 , where (r 2 ) is the mean square atomic radius. (Langevin) 

1 Atoms with a permanent magnetic moment p, have a paramagnetic suscepti¬ 
bility x ~ Np 2 J3'k B T y for pB < k&I, (Gurie-Langevin) 

* For a system of spins S = 4, the exact magnetization is M — 
Nji. tanh(iwhere p — |g(Brilloum) 

* The ground state of electrons in the same shell have the maximum value of S 
allowed by the Pauli principle and the maximum L consistent with this S. The 
J value is L + '$ if the shell is mom than halffell and |L — S | if the shell is less 
than half full, 

* A cooling process operates by demagnetization of a paramagnetic salt at con¬ 
stant entropy. The final temperature reached is of the order of (B JB)T initia j, 
where ife is the effective local field and B is the initial applied magnetic field. 

* The- paramagnetic susceptibility of a '.Fermi ga B- of conduction electrons is 
X =• : .3Np 2 M€p, t: independent of temperature for kgT < e F . (Pauli) 


Problems 

L Diamagnetic susceptibility of atomic hydrogen. The wave function of the hydrogen 
atom in its ground state (Is) is ^ — (7r^|)~ 1/2 exp(—r/oo)> where « 0 = h 2 lme 2 = 
0.529 X T0“ 8 cm. The charge density is p(x t y, z ) = —e|$ 2 according to the statisti¬ 
cal interpretation of the wave function. Show that for this state (r 2 ) = 3a 2 , 
and calculate the molar diamagnetic susceptibility of atomic hydrogen 
(-2.36 X 10" 6 cm 3 / mole). 
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Figure 11 Temperature dependence of the magnetic susceptibility of metals. (Courtesy of C. J. 

Kriessman.) 

2. Hund rules. Apply the Hund rules to find the ground state (the basic level in the 
notation ofTable 1) of (a) Eu + + , in the configuration 4f 5 s 2 p 6 ; (b) Yb 3+ ; (c) Tb 3+ . The 
results for (b) and (c) are in Table 1, but you should give the separate steps in applying 
the rules. 

3. Triplet excited states. Some organic molecules have a triplet (S = 1) excited state at 
an energy J^A above a singlet (S = 0) ground state, (a) Find an expression for the 
magnetic moment (/z) in a field B. (b) Show that the susceptibility for T > A is 
approximately independent of A. (c) With the help of a diagram of energy levels 
versus field and a rough sketch of entropy versus field, explain how this system might 
be cooled by adiabatic magnetization (not demagnetization). 






























4, Heat capacity from internal degrees of freedom* (a) Consider a two-level system with 
an energy splitting k B A between upper-and lower states; the splitting may arise from 
a magnetic field or in other ways, Show that the heat capacity per system is 

„ (W\ , mfe UT 

c “ * (1 + ^7 ■ 

The function is-plotted in Fig, IE. Peaks of this type in the heat capacity are often 
known as Schottky anomalies. The maximum heat capacity is quite high, but for 
T < A and for T P* A the heat capacity Is low, (b) Show that for T ^ A we have 
C = fi; B .(A/iT) 2 + - * * . The hypeifine interaction between nuclear and electronic 
magnetic moments in paramagnetic salts (and in systems having electron spin order) 
causes splittings with A ** I to 100 mK. These splittings are often detected experi¬ 
mentally by the presence of a term in I IT 2 in the heat capacity in the region T ^ A. 
Nuclear electric quadrupole interactions (see Chapter 16) with crystal fields also 
cause splittings, as in Fig, 13. 


I* Pauli spin susceptibility* The spirt susceptibility of a conduction electron gas at abso¬ 
lute zero may be discussed by another method. Let 

N + « *M(1 + 0 ; = iN(l - 0 

be the concentrations of spin~up and spin-down electrons, (a) Show that in a magnetic 
field B the total energy of the spin-up band in a free electron gas is 

e + = e 0 (i + ci 5,3 - mm (i + o , 

where E 0 = in terms of the Fermi energy % in zero magnetic field. Find a 

similar expression for E~. (b) Minimize E to tai — E* + E~~ with respect to f and solve 
for the equilibrium value of f in the approximation £ ^ 1. Go on to show that the 
magnetization if M = 3ZV/tx 2 B/2e F> in agreement with Eq. (42). 

6. €ondmtmn electron ferrmmgmtmm* We approximate the effect of exchange inter¬ 
actions among the conduction electrons if we assume that electrons with parallel 
spins interact with each other with energy —and V is positive, while electrons with 
antiparallel spins do not interact with each other:, (a) .Show with the help of Problem 5 
that the total energy of the spin-up band is 

E + = E 0 (l + 0 sa - £VN 2 (1 + |f - INixB (1 + Cl . 

find a similar expression for E~. (bj) Minimize, the total energy and solve for fin the 
limit l ^ 1. Show that the magnetization is 


M = 


ml B 

2 e F - iVN 


so that the exchange interaction enhances the susceptibility, (c) Show that with B = 0 
the total energy is unstable at f — 0 when V > 4e F i$N* M this is satisfied a-ferromag¬ 
netic state (f ^ 0) will have a lower energy than the paramagnetic state. Because of 
the assumption f ^ 1, this is a sufficient condition for ferromagnetism, but it may 
not be a necessary condition. It is known as the Stoner condition. 
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Figure 12 Heat capacity of a two-level system as a function of T/A, where A is the level splitting. 
The Schottky anomaly is a very useful tool for determining energy level splittings of ions in rare- 
earth and transition-group metals, compounds, and alloys,. 



Figure 13 The normal-state heat capacity of gallium at T < 0.21 K, The nuclear quadrupole 
(C a X" 2 ) and conduction electron (C X) contributions dominate the heat capacity at very low 
temperatures. (After N, E- Phillips.) 
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7, Two-level system. The result of Problem 4 is often seen in another form, (a) If the two 
energy levels are at 4 and “4, show that the energy and heat capacity are 

U - -A tmh{AIk B T) ; C - k B (&/k B Tf seeh 2 (Afk B T) . 

(b) If the system has a random composition such that all values of 4 are equally likely 
up to some limit 4^ show that the heat capacity is linearly proportional to the tem¬ 
perature, provided k & T 4 0 . This result was applied to the heat capacity of dilute 
magnetic alloys by W,. Marshall, Phys, lev,'118, 1519 (1960). It is also used in the 
theory of glasses. Chapter 17. 

8, Paramagnetism ofS = 1 system, (a) Find the magnetization as a function of magnetic 
field and temperature for a system of spins with S = 1, moment /a, and concentration 
n, (b) Show that in the limit /aB ^ kT the result is M = (2npi 2 f3kT)B. 
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